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Abstract 

The D*~ai branching fraction has been measured with data collected by the BABAR 

experiment in 1999 and 2000 corresponding to a total integrated luminosity of 20.6 fb~^. Signal 
events have been selected using a partial reconstruction technique, in which only the af and the 
slow pion (vr^) from the D*~ decay are identified. A signal yield of 18400 =b 1200 events has been 
found, corresponding to a preliminary branching fraction of (1.20 =b 0.07(stat) =b0.14(syst))%. 
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1 Introduction 



A partial reconstruction technique has been used in the past to select large samples of reconstructed 
B mesons with a D*~ in the final state ^ and to measure properties of the B'^ meson. In this 
method the D is not reconstructed, but its four-momentum is inferred from the kinematics of 
the af, the slow pion (vr^) from D*~ decay and the decay constraints. The measurement of the 
branching fraction for B^ D*^a^ is performed as a first step in demonstrating that the partial 
reconstruction method may offer a means of determining the combination of Cabibbo-Kobayashi- 
Maskawa unitarity triangle angles sin(2/3 + 7) using this channej^. 

2 The BaBAR detector 

The BABAR experiment is located at the PEP-II storage ring at the Stanford Linear Accelerator 
Center. A detailed description of the detector and of the algorithms used for the track reconstruction 
and selection of BB events can be found in Ref. For the partial reconstruction analysis of 
B'^ D*^ al only charged tracks are used: particles with transverse momentum pT > 170 MeV/c 
are reconstructed by matching hits in the Silicon Vertex Tracker (SVT) with track segments in the 
Drift Chamber (DCH). Tracks with lower pT do not penetrate a significant distance in the DCH 
and are reconstructed using only the information from the SVT. 

Electron, muon and kaon identification is used in the analysis as a veto in the selection of pions 
forming the ai candidates. Electron candidates are identified by the ratio of the energy deposited 
in the electromagnetic calorimeter (EMC) to the track momentum {E/p) and by the energy loss in 
the DCH (dE/dx). Muons are primarily identified by the measured number of hadronic interaction 
lengths traversed from the outside radius of the DCH through the iron of the instrumented flux 
return (IFR). Kaons are distinguished from pions and protons on the basis of dE/dx in the SVT 
and DCH, and the number of Cherenkov photons and the Cherenkov angle in the Detector for 
Internally Reflected Cherenkov radiation (DIRC). 

3 Data sample 

The data used in this analysis were collected with the BABAR detector in 1999 and 2000. These 
data correspond to an integrated luminosity of 20.6 fb~^ recorded at the T{4S) resonance and, for 
background studies, 2.6 fb^^ collected at 40 MeV below the resonance ("off-resonance" sample). 
Monte Carlo samples of BB and continuum events were reconstructed and analyzed using the same 
procedure as the data. The equivalent luminosity of the generic simulated data is approximately 
one fourth of the on-resonance data, while the number of signal B^ D*~af Monte Carlo (MC) 
events is approximately ten times the number expected in the on-resonance data. 

4 The partial reconstruction technique 

In the partial reconstruction of the decay chain B^ D*^af, followed by D*~ D^tt^, only 
the ai and the vr^ from D* decay are required. For this analysis the ai is reconstructed via the 

'^Since the selection of D*^a^ and D*'^a1 are identical, the charge conjugate state is implied throughout the 
paper. 
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decay chain af — > p^ir^. The angle between the momentum vectors of the B and the oi in the 
center-of-mass frame (CM) is then computed: 

, Ml,- - M|o - Ml + EcMEg, 

cos OBai = np \- \ (1) 

where is the mass of particle x, Ea-^ and are the measured CM energy and momentum of 
the ai, EcM is the total CM energy of the beams and Pb = \J E'^j^/A — M^o • Given cos^sai and 
the measured four-momenta of the tTs and the oi, the B four-momentum can be calculated up to 
an unknown azimuthal angle (j) around ■ For every value of (p, the expected D four-momentum, 
P-j^o, is determined from four-momentum conservation and the (/)-dependent "missing mass" is 

calculated, m^iss{4>) = ^\Pj^o p. Defining nimin and nimax to be the minimum and maximum 
values of mmiss(</') obtained by varying (f) over the allowed range, the missing mass is calculated 
as m-miss = ^(mmax + mmin)- For signal events, this variable peaks at the mass, while for 
background events it has a broader distribution. For this reason, mmiss can be used to determine 
the fractions of signal and background events in the data sample. 

5 Event selection 

Data and Monte Carlo events are first selected with the following loose requirements: 

• R2, the ratio of the 2nd to the 0th Fox Wolfram moments Q, less than 0.35; 

• at least one ai candidate^ such that: 

- the oi invariant mass nia^ is between 1.0 and 1.6 GeV/c^; 

- the oi momentum pa^, computed in the CM frame, is between 1.85 and 2.30 GeV/c; 

- the vertex probability obtained from a vertex fit of the 3 pions is greater than 1%; 

- the invariant mass of at least one of the two possible 7r"^7r~ combinations (w-^tt) is in the 
range [0.278, 1.122] GeV/c^; 

• at least one additional track (the vr^) with CM momentum p^^^ between 50 and 700 MeV/c; 

• for each selected D*ai candidate, there must be at least 2 additional particles (charged or 
neutral) . 

The fraction of events selected by applying these cuts is summarized in Table |^ for signal and 
background MC samples and for off-resonance data. 

The main source of background in this analysis is continuum qq events, where q = u,d, s,c. A 
neural network is used to optimize the separation of BB events from the continuum background, 
independent of any particular B decay channel. The NN has three layers, with 11 input nodes, 
15 hidden nodes and one output. Its definition relies on the different topologies of signal and 
background at the T{4S): while BB events are more "spherical", qq events are more "jet-like". 
The variables used to discriminate jet-like from isotropic events are R2; the thrust of the event Q; 
the two invariant masses squared, obtained by adding the four momenta of all particles going into 

^By we refer to the tt"'" tv~ tt"*" final state. 
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each of the hemispheres divided by the plane perpendicular to the thrust axis; and the angle of the 
thrust axis with respect to the e~^e~ direction. 

The signal to background separation becomes more complicated when one or more gluons are 
emitted, which is more likely to happen when light quarks are produced. In this case there are 
multiple preferred axes, so that the event shape more closely resembles that of the signal. To 
discriminate further between signal and background in this case, the tracks are clustered into 3 or 
4 jets using the so called "Durham" algorithm [^,^. 

The event is first clustered into four jets and the following discriminating variables are added 
to those previously defined in order to build up the network: y4, the jet metric [^,^ obtained when 
the event is clustered in 4 jets; the QCD matrix element |^; the cosine of the maximum angle 
between each pair of jets; the angle between the plane defined by the two most energetic jets and 
that defined by the other two jets; and the angle between the lowest and second lowest energy jet. 
The event is then further clustered into 3 jets and the jet metric ys is also added to the list of 
network input variables. 

The distributions of BB MC events and on-resonance data that pass the event selection criteria 
are shown in Fig. |^. The on-resonance distribution is shown after the subtraction of the off- 
resonance distribution, scaled by the ratio of the on-resonance to off-resonance luminosities and 
the CM energy squared. By selecting events for which the NN output (Onn) is greater than 0.25, 
66% of the continuum events are rejected. The fractions of events selected by this requirement in 
the various analyzed samples are shown in Table |l[ 

After the NN cut, further selection requirements are applied. There must be at least one oitt 
combination with net charge equal 0; the three pion invariant mass must be in the range [1.0, 
1.45] GeV/c^; the three pion momentum in the CM frame must be in the range [1.9, 2.25] GeV/c; 
at least one 7r"^7r~ combination from the three daughters of the oi must have invariant mass in the 
range [0.65, 0.9] GeV/c^; and the vertex probability of the 3tt plus the tTs must be greater than 3%. 

The fractions of events selected by these requirements are summarized in Table ||. 

Table 1: Fractions of events in MC data that are selected by the requirements applied in the 
analysis. 



Cuts applied 


B° D*-ai 


B^B^ 


B+B- 


qq 








(no signal) 




MC 


Reconstruction Efficiency 


38.: 


i% 


14.6% 


13.9% 


6.4% 


Onn > 0.25 


85.4% 


91.9% 


92.1% 


45.2% 


Qai-TT = 


94/. 


i% 


86.7% 


85.2% 


83.5% 


1.0 < < 1.45 GeV/c^ 


81.7% 


71.8% 


71.7% 


73.6% 


1.9 < Pa, < 2.25 GeV/c 


94.^ 


i% 


92.4% 


92.5% 


93.3% 


0.65 < m^^ < 0.9 GeV/c2 


76.^ 


1% 


59.1% 


60.5% 


64.3% 


Vertex Prob in > 0.03 


80.( 


)% 


74.2% 


73.9% 


76.9% 


Total Efficiency 


14.( 


3% 


3.4% 


3.2% 


0.82% 
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Figure 1: Distribution of the neural net output for (a) B'^B^ MC, (b) B^B^ MC and (c) on- 
resonance (off-resonance subtracted) data events, with superimposed distribution for all BB MC 
events. All events satisfy the event selection criteria. 
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6 Results 



Applying all of the selection criteria, the ?n.miss distribution of on-resonance and off-resonance data 
is obtained for "right-sign" events, which are events where the ai and tTs candidates have opposite 
electrical charges. The off-resonance distribution is then scaled to take into account the difference 
in luminosities and CM energies between the two samples and subtracted bin-by-bin from the on- 
resonance distribution. The resulting plot is fitted, using a minimum technique, to a linear 
combination of the m^[ss distributions of: 

1. BB Monte Carlo events, excluding correctly reconstructed signal events; and 

2. correctly reconstructed signal Monte Carlo events. 
In Fig. 1(a) the 

^miss distribution of "right-sign" on-resonance data, off-resonance subtracted, 
is shown, together with the distributions of BB background MC and signal MC events. The signal 
yield that is obtained is 18400 ± 1200 events. The BB contribution to the fit is 0.995 ± 0.015 
of the value expected, given the total number of BB events in the data and in the Monte Carlo 
simulation. 

Given this yield, the oi branching fraction]^, the total signal efficiency of Table Q and the fraction 
of events with multiple signal candidates, the following preliminary result for the B^ —>■ D*^ai 
branching fraction is obtained: 

B{B^ D*-a+) = (1.20 ± 0.07)%, (2) 

where the error is statistical only. The result is in very good agreement with the current best 
measurement of (1.30 ± 0.27)% ||-|n|. 



Several tests were conducted to verify that the background shapes in data and MC agree 
and do not give rise to a spurious signal. Charged ai candidates were combined with tracks of 
the same charge into "wrong-sign" combinations, and were analyzed in the same way as "right- 
sign" combinations. The fyimiss 

distribution of these candidates, shown in Fig. |2|(b), shows a good 
agreement between data and MC in the signal region {m^iss > 1.854 GeV/c^). The BB contribution 
to the fit is 0.990 ± 0.016 of the value expected from MC simulation. 

Among events collected on-resonance, the fraction containing more than one aivr combination 
that passes all the analysis requirements in the signal region, defined by rriYniss > 1-854 GeV/c^, 
is F = (15.10 ± 0.14)%. This agrees with the fraction Fmc = (15.51 ± 0.26)% obtained from a 
weighted mix of off-resonance, BB and signal MC events. It was verified that the reconstruction 
efficiency in qq MC and off-resonance data are in good agreement. 

7 Systematic errors 

The selection criteria applied in the analysis are varied within a reasonable range around their 
chosen values. The branching fraction and its error are recalculated for each value, obtaining N 



different measurements for N different choices of the requirement [12|. The variation of the N 



results with respect to their average, taking into account statistical correlations between them [13|, 
is taken as a systematic error. The method of Ref. [14| is used in order to disentangle the statistical 
fluctuations from the systematic ones. 

^In this analysis it is assumed that B{a^ — + p"^^) = 0.4920, based on isospin considerations and phase space 
corrections. 
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Figure 2: (a) mmiss distribution of continuum-subtracted on-resonance data events (data points), 
BB background MC events (dashed histogram) and BB background plus signal events (solid his- 
togram) for "right-sign" oivr combinations. The histograms are the result of the fit procedure 
described in the text, (b) Same distributions for "wrong-sign" oivr combinations. 
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The final contributions to the error on the branching ratio due to the variation of the selection 
requirements are shown in Table |2|. It has been verified that the systematic error obtained for each 
requirement does not depend on the variation studied. 



Table 2: Systematic errors on B{B^ 
ments in the analysis. 



D* 



due to changing the value of some of the require- 



Cut applied 


Range 


Studied 




N 


Error (%) 


OnN > Omin 


Oniin ~ 


: 0.15 to 0.35 


1.37 


5 


4.2 




^min ~ 


= 1.0 GeV/c2 










^max " 


= 1.4 to 1.5 GeV/c2 


0.99 


5 


negl. 


I'min ^ Pai ^ Pmax 


Pmin — 


1.85 to 1.95 GeV/c^ 










Pmax — 


2.15 to 2.25 GeV/c^ 


1.96 


5 


6.6 


Vertex Prob Svr > Pmin 


P . — 

^ mm 


0.01 to 0.15 


1.15 


5 


2.4 


Vertex Prob Att > Pmin 


P ■ — 

^ mm — 


0.01 to 0.1 


1.10 


5 


2.2 




^min ~ 


= 0.65 to 0.75 GeV/c^ 










^max " 


= 0.8 to 0.9GeV/c2 


1.45 


5 


4.6 


Total error 








9.6 



A conservative systematic error of 0.35% is determined from MC due to the dependence of the 
TTs reconstruction efficiency on the tTs momentum. The systematic error due to track reconstruction 
efficiency is 4.2% and the uncertainty in the total number of B mesons in the data sample is 1.6%. 

The BE background has a component that peaks slightly under the signal region (Fig. |2|). This 
is due mostly to signal events in which one or more of the selected tracks did not originate from the 
signal B^. The contribution of this background is varied in the BB MC by ib-\/(0.07/1.2)-^ + 0.104-^, 
i.e., by the relative statistical error in the central value of the branching fraction plus the total 
relative systematic error calculated up to this point, added in quadrature. This results in a 4.5% 
variation of the signal yield, which is added to the total systematic error. The total systematic 
error is 11.5%. 

Other possible sources of systematic error have been investigated. The non-resonant decay 
channel 5° D* 

ing fraction for this mode is 0.57 ± 0.31% 



~ p^TT~^ could contribute to the quoted branching fraction. The measured branch- 



10 1 . Since the central value is inconsistent with the 



total B{B^ D* ir^TT ir~^), this channel is ignored in our fit. As a result, there is a potential 
contribution from this channel that is included in the quoted branching fraction for B^ — s- D*~at . 



D* 



at a rate of 



The central value will shift according to the branching fraction for B^ 
-3.3% X ^ D*- A+)/0.57%. 

Likewise, the decay B — > D**ai could affect the signal yieldQ. To study its effect, D** MC 
events are added to the generic BB sample at the level of B** = B{B D**ai) x B{D** — > 
Z)*7r) = 0.35% [|^, and the fit to the missing mass distribution is repeated. This results in a 
reduction of the signal yield of 4.3%. Since the B{B — > D**ai) has not yet been measured, based 
on this result the branching fraction obtained in this analysis will be shifted by contributions from 
B D**ai at a rate of -4.3% x e**/0.35%. 



'D** denotes the sum of Di, and D2 states. 



14 



8 Conclusions 



With a partial reconstruction technique, 18400 ± 1200 i? D* af events have been found in the 
BABAR data set of 20.6 fb~^ on-resonance events. This corresponds to the branching fraction 

B{B° D*-a^) = (1.20 ± 0.07 ± 0.14)%, (3) 

where the first error is statistical and the second is systematic. Due to the uncertainty in the contri- 
bution of B^ —> D*^p^TT^ events in the signal sample, the central value of the quoted branching frac- 
tion will shift according to the branching fraction for B^ — > D*' at a rate of —3.3% x B{B 
L>*-p°7r+)/0.57%. Likewise, due to the unknown value of B** = B{B D**ai) x B{D** D*7r) 
the central value of the branching fraction will be shifted by contributions from B — > D**ai events 
at a rate of —4.3% x i3**/0.35%. The result is in good agreement with the current world average 
value of (1.30 it 0.27)% [P-[Tl| but reduces the uncertainty by a factor of two. 
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